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Abstract. High resolution spectra with UVES/VLT of SN 1987A from December 2000 until November 2005 show 
a number of high ionization lines from gas with velocities of ~ ±350 km s~^, emerging from the shocked gas 
formed by the ejecta - ring collision. These include coronal lines from [Fe X], [Fe XI] and [Fe XIV] which have 
increased by a factor of ~ 20 during the observed period. The evolution of the lines is similar to that of the soft 
X-rays, indicating that they arise in the same component. The line ratios are consistent with those expected from 
radiative shocks with velocity 310 — 390 km s~^, corresponding to a shock temperature of (1.6 — 2.5) x 10^ K. A 
fraction of the coronal emission may, however, originate in higher velocity adiabatic shocks. 
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1. Introduction 

The collision of the ejecta of SN 1987A with the 
circumstellar ring has been obse rved in many wave- 
length ranges (e.g.. iMcCravl l)2005|) for a review). In the 
UV/optical range this has been monitored with HST 
showing an increasing nuinber of hot spots ar ound the ring 
l)Michael et al.ll2000ll20o3 |Pun et al.ll2002tl . Most likely, 
these hot spots are caused by the impact of the blast wave 
on protrusions from the inner circumstellar ring. As time 
progresses an increasing number of these are expected to 
appear, and finally the whole ring will be immersed in the 
collision. 

The spectroscopic HST observati ons have identified 
three different velocity components (P un et alJ2 002'). One 
narrow, nearly Gaussian, with FWHM 10 km s~^, com- 
ing from the unshocked ring gas. Secondly, there is an in- 
termediate velocity component arising from shocked cool- 
ing gas, with a clearly non-Gaussian form, extending to 
~ ±300 km s~^. Finally there is a very broad component 
extending to 15 000 km s""'^, probably coming from the 
reverse shock, resu lting from the interaction with the sur - 
rounding medium l)Michael et al.ll998tlSmith et al.l20 0^. 
We will refer to these as the narrow, intermediate and high 
velocity components, respectively. 



The collision is also seen in the radio as a ris- 
ing non-thermal svnchrotron flux l|Gaensler et alJ Il997t 
iManchest er et all l200i |2005|) . The fight curves in this 
range sh ow a roughl y linear increase with time. Further, 
iBouchet et al.l ()2006|) find from ground based and Spitzer 
mid-IR observations evidence for hot dust as well as highly 
ionized line emission. 

In X-rays ROSAT and Chandra have observed SN 
1987A at several epochs, and also in this energy range a 
steadily rising flux is seen jPark et al.ll2004ll2005|) . While 
the hard X-rays correlate well with the radio flux, the rise 
time in the soft X-rays is considerably shorter. Of spe- 
cial i nterest are the high resolution spectral observations 
(Zhc kov et aLll2005lEl0l . which show a number of lines 
from H and He like ions of N, O, N e, Si, Mg and S, as well 
as Fe XVII. IZhekov et all (l2005l) argue that these lines 
arise in the radiative shocks propagating into the protru- 
sions, as well as the reflected shock propagating back from 
these. In this paper we discuss optical observations com- 
plementary to these X-ray observations of highly ionized 
gas. 

In iGroningsson et ahl 



Send offprint requests to: P. Groningsson 

* Based on observations performed at 
Southern Observatory, Paranal, Chile. 



the European 



l|200()|) we report on spec- 
tral shapes and fluxes of optical lines created by the 
ejecta/ring interaction. As we use high-resolution echelle 
spectroscopy, we can easily disentangle the different ve- 
locity components from each other. Here we concen- 
trate on the evolution of intermediate velocity high- 
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ionization lines, such as [Fe X] A6374.5, [Fe XI] A7891.8 
and [Fe XIV] A5302.9, as they are efficient probes of the 
shocked gas. 

In our analysis we use a recession velocity of 
281 km s~^ (the bulk motio n of the northern side o f 
the inner circumstellar ring ((Gro ningsson et al. "2006^) 
and a reddening of Eb-v = 0.16 |Fitzpatrick & Walborn 
19901) with Eb-v = 0.0 6 from the Milky way (e.g., 
Stavelev-Smith et alJl2003|) and Eb-v = 0.10 from the 
LMC. The reddening law was taken from ICardcUi et al] 
l|l989|) using i?y = 3.1. The differences between the LMC 
extinction law and the Galactic law are negligible in the 
optical at low color excess l|Fitzpatrick 1999) and have 
therefore been ignored. In Sect. 2 and 3 we describe the ob- 
servations and observational results, respectively. In Sect. 
4 we discuss modeling of the lines, followed by a general 
discussion in Sect. 5 and conclusions in Sect. 6. 

2. Observations 

SN 1987A was observed in service mode with the 
Ultraviolet and Visual Echelle Spectrograph (UVES) 
at ESO/VLT at Paranal, Chile. UVES is a cross- 
dispersed echelle spectrograph covering (with two differ- 
ent dichroic settings ) the spectral range ~ 3100-10000 A 
i Dekker et al The light beam is split up in two sep- 

arate arms. The arm covering the blue part of the spec- 
trum (~ 3100 - 4900 A) has a single CCD detector with 
a spatial resolution of O'.'246/pix. The other arm cover- 
ing the red part (~ 4800 - 10000 A) is equipped with 
two CCDs having the somewhat higher spatial resolution 
of O'.'182/pix. The slit width used resulted in a resolving 
power of ~ 50 000 corresponding to a spectral resolution 
of ~ 6 km s^^. 

Spectra of the ring and ejecta of SN 1987A were ob- 
tained on October 16 1999, December 10-14 2000, October 
4-7 2002, March 21 - April 12 2005 and October 20 - 
November 12 2005, in the following referred to as epochs 
1-5, respectively. In all cases, except for epoch 1, the po- 
sition angle was PA= 30°. For epoch 1 it was PA= 20°. 
The log of the observations is given in Table The ob- 
servations and data reduction are discussed in detail in 
iGroningsson et a~ ] §006). As can be seen from Table ^ 
the exposure of the first epoch was very short and the S/N 
low. This made it impossible to obtain meaningful fluxes 
for the coronal lines for this epoch, and therefore we do 
not discuss it further in this paper. 

Because the spatial resolution of these ground-based 
observations is limited compared to the dimensions of the 
ring, we cannot distinguish between different hot spots lo- 
cated on the same side of the ring. To identify how many 
hot spots were covered by the slit we retrieved and stud- 
ied HST images taken with the WFPC2 and ACS instru- 
ments at roughly the same epochs. These show that only 
Spot 1 is responsible for the emission from the shocked 
ring at epoch 1. At epoch 2 Spot 1 still dominates the 
shocked gas emission on the north side of the ring, but 
also two spots on the opposite side of the ring are now 



prominent. At epoch 3, three to four different hot spots 
contribute to the emission on the northern side and by 
the time of the last epochs the entire ring is lit up by the 
ejecta-ring collision. Because of the difference in velocity 
and to isolate a limited part of the ring we will in the fol- 
lowing only discuss the spectrum from the north side of 
the ring, i.e., where Spot 1 first appeared. A comparison of 
th e the kinematics and flux of the two sides will be given 
in lCroningsson et al. 1 1|2006|) . 

To estimate the accuracy of the absolute flux level of 
the spectra we compared the individual sensitivity func- 
tions created from the reduced standard star spectra and 
their corresponding physical fluxes. The sensitivity curves 
differed typically by ^ 10% from each other. These ex- 
posures of the standard stars used a wide slit (10") and 
hence no slit losses are likely to occur. Our science data 
exposures, on the other hand, used a relatively narrow slit 
and since SN 1987A is an extended source the slit losses 
could be considerable and depend on atmospheric condi- 
tions such as the seeing. To estimate how much the atmo- 
spheric conditions could influence our results we have com- 
pared line fluxes of the individual science spectra for the 
epochs which have more than one exposure. This compar- 
ison reveals that the line fluxes differ by ^ 10% from one 
exposure to another. Finally, to estimate the total system- 
atic error, we have co mpared our fluxes with H ST spectra 
and photometry C see IGroningsson et al. lEoO^ . and find 
that the accuracy of the VLT fluxes should in general be 
better than - 30%. 

3. Results 

In IGroningsson et ahl l)2006l) we give a full list of the 
different lines present, both from the unshocked and the 
shocked ring. The former are dominated by thermal broad- 
ening with a FWHM of ~ 10 km s~^. The lines from 
the shocked gas are on the other hand dominated by the 
macroscopic velocity and are highly asymmetrical, extend- 
ing to ~ 250 - 400 km s~^ (HWZI). The two components 
can therefore easily be well separated, and a subtraction 
of the broader from the narrower component by interpola- 
tion of the intermediate components by regression splines 
is fairly straightforward. More problematic is the blending 
of different lines from the shocked component. In partic- 
ular, the large number of Fe II lines from this component 
causes a considerable number of blendings. 

The fuU spectrum between 3100-10000 A contains 
~ 170 intermediate velocity lines. About one third of 
these are Fe II lines, but there are also strong lines from 
more highly ionized ions, like [O III] AA4958. 9, 5006.8, 
[Ne III] AA3868.9, 3967.5, [Ne V] AA3345.3, 3425.9, and 
[Ar III] A7135.8. The most interesting result is, however, 
that we find a number of highly ionized coronal lines, 
which can only be seen thanks to the high S/N and high 
spectral resolution of our UVES observations. 

In Figs, n - 01 we show a compilation of the 
most important of these lines for the December 2000, 
October 2002, March/ April 2005 and November 2005 
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Table 1. Log of UVES observations. 



Epoch 


Date 


Days after 
explosion 


Grating 


A range 

A 


Slit width 
arcsec 


Resolution 
A/AA 


Exposure 
s 


Seeing 

arcsec 


1 


1999 Oct 16 


4618 


DICl 346+580 


3100-3900 


1.0 


40,000 


1,200 


1.0 










4800-6800 










2 


2000 Dec 9-14 


5039-5043 


DICl 346+580 


3100-3900 


0.8 


50,000 


10,200 


0.4-0.8 










4800-6800 














5038-5039 


DIC2 390+860 


3300-4400 






9,360 


0.4 










6700-10000 










3 


2002 Oct 4-7 


5704-5705 


DICl 346+580 


3100-3900 


0.8 


50,000 


10,200 


0.7-1.0 










4800-6800 














5702-5703 


DIC2 437+860 


3800-4900 






9,360 


0.4-1.1 










6700-10000 










4 


2005 Mar 21 


6601-6623 


DICl 346+580 


3100-3900 


0.8 


50,000 


9,200 


0.6-0.9 




Apr 8-12 






4800-6800 














6621-6623 


DIC2 437+860 


3800-4900 






4,600 


0.5 










6700-10000 










5 


2005 Oct 20 


6826 


DICl 346+580 


3100-3900 


0.8 


50,000 


2,300 


0.9 




Nov 1-12 






4800-6800 














6814-6837 


DIC2 437+860 


3800-4900 






9,200 


0.5-1.0 



6700-10000 



epochs. These include [Ne V] A3425.9, [Ar V] A7005.7, 
[Fe VII] A6087.0, [Fe X] A6374.5, [Fe XI] A7891.8, and 
[Fe XIV] A5302.9, which are all likely to originate 
from the cooling post-shock gas. In addition, we include 
for comparison the [O III] A5006.8, [Ne III] A3868.8, 
[Fe II] A7155.2, [Fe III] A4658.0 lines, which are, however, 
likely to come from considerably cooler, photoionized gas 
(Sect. 01). Because these are very useful as diagnostics of 
the shocked ring g as we concentrate in t his pa per on these 
lines, and refer to iGroningsson et al. I l)2006(l for a more 
detailed analysis of the full spectrum. 

As an important example we show in Fig. El the re- 
gion of the [Fe XIV] A5302.9 line. It is here seen that the 
blue wing of the [Fe XIV] A5302.9 line is blended with the 
intermediate velocity [Fe II] A5296.8 line. To remove this 
we use the [Fe II] A5527.6 line as a template. The rea- 
son why we have chosen this line instead of the stronger 
[Fe II] A7155.2 line is the similar excitation potential to 
the [Fe II] A5296.8 line. The result of this subtraction can 
be seen in the lower panel of Fig. |31 In addition, we also 
see a narrow component of [Ca V] A5309.1 at 350 km s~^. 

The [Fe X] A6374.5 line is close to the red wing of the 
very strong [O I] A6363.8 line. In addition, there is a nar- 
row Si II A6371.4 hne close to its peak. The [O I] A6300.3 
line was used as a template to deblend [O I] A6363.8 from 
the [Fe X] A6374.5 line. The extension of the blue wing is 
therefore for this line uncertain. 

In Fig. El we show the most important high ionization 
lines from the October 2002 observations after deblend- 
ing. The other epochs have been deblended in the same 
manner. 



We have also searched the spectra for [S XII] A7611.0, 
[Ar XIV] A4412.2 and [Ca XV] AA5446.4, 5694.4. The 
[S XII] line unfortunately falls on top of some very strong 
atmospheric bands. While in the data there is some evi- 
dence that a line may be present at the wavelength of the 
[S XII] line we cautiously attribute this to the effects of 
the atmosphere. We have also searched HST spectra in 
this spectral region. Because of the low dispersion, 4.92 A 
pixel"^ or 194 km s~^ pixel"^, and fairly low S/N, only a 
weak upper limit of ^ 10^^^ erg s~^ for the flux around 
our epoch 3 can be set. 

In the region of the [Ar XIV] A4412.2 line there is a line 
redshifted from this wavelength by ^ 200 — 300 km s^^ 
in our VLT spectra. This is, however, likely to be a blend 
of [Fe II] AA4413.8, 4416.3, which matches weh with the 
velocity of the observed feature. Subtracting these lines, 
using [Fe II] A7155.2 as a template in the best S/N spec- 
trum of this region at epoch 5, we estimate an upper 
limit of ^ 15% of the [Fc XIV] flux at this epoch, i.e. 
< 2.2 X 10-1^ erg s~icm"2. 

We have also looked for the [Ca XV] A5694.4 hne, 
but do not detect this to a flux limit of ^ 15% of 
[Fe XIV] A5302.9 at epoch 4. 

In Table (21 we give the fluxes of the intermediate ve- 
locity component for the highly ionized lines in the spec- 
trum at the four well observed epochs. These fluxes are 
corrected for blending, as described above, but not for 
reddening. The estimated extinction correction factors 
are listed in table El In addition, we detect the [Fe VI] 
A5335.2 line at epoch 4 with a flux of (1.3 ± 0.2) x 
10"^® erg s^^cni^-^. 
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Fig. 1. Original line profiles from December 2000 (our 
epoch 2) of the most important high ionization lines from 
the north side of the shocked ring. The [Fe III] A4658.1 
line was not covered by the observations at this epoch. 
Note the different flux scales. All velocities are with re- 
spect to the systemic velocity of the northern part of the 
ring, 281 km s~^ 



As we see from Table 12 and Figs. ^ - 0] there is an 
increase in the fluxes of all lines between epochs 2 and 5. 
Over this period the [Fe X], [Fe XI], and [Fe XIV] lines 
have increased by a factor ~ 24 — 30. At epoch 2 some 
of the weaker lines are, however, close to the noise level, 
and only the fluxes at epoch 3 are reliable. In Fig. [7| we 
show the evolution of the fluxes at the observed epochs. 
To compare with the X-ray and radio flux evolution, the 
fluxes in the flgure have been normalized to the October 
2002 level, i.e., epoch 3. From this Figure it is seen that the 
[Fe X] , [Fe XI] and [Fe XIV] lines evolve very similar to the 
flux evolution of the soft X-rays, and considerably faster 



Fig. 2. Same as Fig.^ but for October 2002. 

than the hard X-rays or the radio. The [Ne V], [Ar V], and 
[Fe VII] lines, however, increase considerably more slowly. 
We return to this result in Sect.El 

Although there seem to be considerable differences be- 
tween the line proflles of the different ions, it is impor- 
tant to realize that both the deblending procedure and 
the low S/N for the faintest lines may introduce con- 
siderable uncertainties. Only for lines of high S/N is it 
meaningful to compare the line profiles directly. The ve- 
locities of the different lin es will be discussed in detail in 
ICroning-sson et all (j5 



Here we only make some brief 
remarks about the most relevant issues. If we compare the 
red extension of [Fe XIV] A5302.9 with the [Fe II] A7155.2 
line we find evidence for an extension to 350 km s~^ for 
[Fe XIV], while the [Fe II] fine extends to - 250 km s^^. 
There are, however, important systematics due to uncer- 
tainties in background subtraction both in the spatial and 
dispersion directions. In addition, deblending of line pro- 
files gives rise to uncertainties especially when the de- 
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Fig. 3. Same as Fig-H but for March - April 2005. 



Fig. 4. Same as Fig.^ but for November 2005. 



blended line is strong compared to the line in consider- 
ation. This makes e.g., the extent of the blue wing of 
the [Fe XIV] A5302.9 line more uncertain than the red. 
With the systematic uncertainties included we estimate 
that the values for maximum velocities should be accu- 
rate to within 50 km s~^. 

4. Origin of the coronal emission 

The intermediate component of both the [O III] A5006.8 
and the [Ne V] A3425.9, as well as several 
strong UV lines, like N V AA1238.8, 1242.8, 
Si IV AA1393.8, 1403.8, C IV AA1548.2, 1550.8, C II A2325 
a nd Mg II AA2795.5, 2802.7, were seen in the HST spectra 
of IPun et al.. (j2002) . However, no coronal lines were seen 
probably due to the lower spectral resolution, as well as 
the lower S/N. The [Fe X] A6374.5 and [Fe XIV] A5302.9 
coronal lines have, however, been observ ed in several su- 
pernova remnants, like the Cygnus loop ijWoodgate et alJ 



Il974l). Pun A. RCW8 fi llT.ucke et a,lJll979ll and IC433 
JWoodga,te et ah 1979). 

IPun et al.l (j2002'l propose that the intermediate com- 
ponent arises when the blast wave hits dense protrusions 
from the ring. The velocity of the shocked gas behind the 
transmitted shock is given by ~ (pHii/pspot)^^^14iast- Here 
Phii is the density of the H II regi on of the progenitor star 
Estimates give pnii ^ 10^ cm~^ l)Chevalier fc DwarkadasI 
Il99.^ . while the density of the hot spot should be close 
to that of the main ring materia l ^ 10'* cm~^. For a blast 
wave velocity of ~ 3800 km s"* jPark et alJl2005|) the ve- 
locity of the transmitted shock should be ~ O.lVbiast ~ 
300 — 500 km s^*. This is close to that observed from the 
intermediate component. 

The origin of the X-ray emission is, however, not so 
clear. Part of the X-rays may arise from the blast wave 
propagating with a velocity of ^ 3 500 km s~* into the 
H II region of the proge nitor, interior to the equatorial 
ring. IZhekov et al.l l)2005|) . however, find a velocity from 
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Table 2. Fluxes (uncorrected for reddening) of the high ionization lines from the north side of the shocked ring. 





Rod" ^ll7■c^■\7dl 




Flux 






A 




lO"^** erg s"^cm"^ 


correction 


fNe VI 

[INC VJ 


3425.86 


2 




2.06 






3 


5.25 ±0.61 








4 


8.58 ± 1.26 








5 


9.69 ± 1.83 




fAr VI 


7005.67 


2 




1.41 






3 


0.34 ±0.07 








4 


0.31 ±0.09 








5 


0.35 ±0.08 




[Fe VII] 


6087.0 


2 




1.50 






3 


0.80 ±0.11 








4 


1.23 ±0.15 








5 


1.50 ±0.26 




\Fe XI 


6374.51 


2 


0.51 ± 0.11 


1.47 






3 


3.11 ±0.24 








4 


6.54 ±0.37 








5 


9.27 ±0.69 




[Fe XI] 


7891.94 


2 


0.91 ±0.17 


1.32 






3 


2.98 ±0.25 








4 


9.10 ±0.40 








5 


9.65 ± 0.34 




[Fe XIV] 


5302.86 


2 


0.67 ±0.12 


1.61 






3 


3.89 ±0.31 








4 


11.00 ±0.39 








5 


14.89 ± 0.69 





the line profiles much smaller than that inferr ed from the 
expan sion from the size o f the X-ray image llPark et 8,1.1 
I2OO4I) and the radio source ((Manchester et al.l2n02tl20n,'il) . 
which implies an expansion rate of ~ 3 000 km s~^. 

Zhekov et al. therefore propose a scenario where the 
blast wave propagates with this comparatively high ve- 
locity inside the ring in a medium which is dominated by 
a density characteristic of the H II region, ^ 10^ cm~'^. 
In this scenario most of the X-ray line emission arises 
where the blast wave hits the dense protrusions from the 
ring, forming tr ansmitted shocks i nto th e dense gas. In 
addition to this, iBorkowski et al. find that there 

will be reflected shocks going back from the hot spot 
into the shocked gas behind the blast wave. The ve- 
locity of the shocked gas behind the reflected shock is 

(pHIl/pblast)'/'Klast « O.^Vue^t « 1 000 - 1 500 km S"! . 

This, and the even lower velocity from the transmitted 
shocks, explains the low velocity compared to the blast 
wave inferred from the X-ray line emitting gas. The higher 
density behind reflected shocks, as well as the transmit- 
ted, causes these to dominate the X-ray emission over the 
blast wave. It also suggests that most of the X-ray emis- 
sion should be correlated with the optical emission from 
the hot spots. Because of the lower density and higher 
temperature behind the reflected shocks, these are likely 
to be adiabatic. The cooling gas seen in the coronal lines, 
as well as in the optical/UV emission from lower ionization 



stages, is therefore likely to arise in the radiative shocks in 
the hot spots. We will therefore investigate this possibility 
more quantitatively below. 

For a gas with x = n(IIe)/n(II) — 0.25 
ijLundavist fc FranssonlUfjQd) the temperature behind the 
shock going into the hot spot is 

9(1 + 42^), 



2.27 X 10" 



1.5 X 10" 



(2- 



-3a;) 



300 km s-i ' ^' ^'^^ 

iMichael et alJ l)2002j) find that for a shock with veloc- 
ity 500 km s"'^ equipartition between electrons and ions 
takes place on a time scale of ^ 60 (rie/lO* cm~'^)~^ days. 
Shocks with lower velocity will equilibrate faster. In the 
following we therefore assume that To = Tjon = Tg. 

For 100 ^ 14 ^ 600 km s~^ the cooling func- 
tion fo r the CNO-enriched c o mposi tion found for the 
ring in iLundavist & Franssoi] l(l996l) is A(To) w 4.0 x 



lO-^^lTjW K)-°-J erg s-^cm^, using the code by 
iNvmark et alJ l)2005f) . The cooling time of the shock may 
therefore be approximated by 



cool 



f-spot 



104 



1-3 



300 km s-i 



3.4 



years. 



(2) 



where ngpot is the pre-shock density. Therefore, depending 
on the density and the shock velocity, the shock may be 
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Fig. 5. Top panel: Line profile of [Fe XIV] A5302.9, 
before debelnding. Middle panel: Template profile for 
[Fe II] A5296.8, based on the [Fe II] A5527.6 line. Bottom 
panel: Deblended profile. The narrow line at 350 km s^^ 
in the top and bottom panels is [Ca V] A5309.1. 

either adiabatic or radiative. While coronal emission lines 
may arise in any hot gas with a temperature of ^ 10® 
K, the wide range of ionization stages with similar line 
profiles seen in the intermediate velocity component, from 
Fe II to Fe XIV, strongly indicates that at least some of 
the emission is originating in radiative shocks. 

The structure of the radiative shocks was discussed in 
|Pun et a l. (2002). Pun et al. , however, mainly focused 
on low and intermediate ionization lines. To calculate the 
flux of the optical/UV/IR coronal lines, as well as some 
of the other high ionization lines, we com pute the spec- 
trum of the shock, using the shock code of iNvmark et alJ 
l)2005|) . Although originally used for the reverse shock in 
supernovae, this is, of course, applicable for any radiative 
shock, independent of shock direction. This determines the 
temperature and ionization structure of the shocked gas 
by solving the time dependent ionization and hydrody- 
namic equations, assuming that the structure is station- 
ary and one-dimensional. The latter should be sufficient as 
long as the shock is radiative. CoUisional ionization rates, 
as well as recombination rates, ar e from the mos t recen t 
data, and references are given in iNvmark et al.l lj200fil) . 
In addition, we include charge transfer, which becomes 
important at these comparatively low temperatures. The 
uncertainties of the ionization and reco mbination dat a es- 
pecially for iro n have been discussed and 
ICianetti et alJ l)200nl) . The emissivity is calculated using 
full multi-level atoms. Pre-ionization of the unshocked gas 
by X-rays from the shock is calculated using an updated 
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Fig. 6. Deblended line profiles from October 2002 of the 
high ionization lines in Fig. [21 



version of the code in IChevalier fc FranssonI ll 19941). Th e 
abundances used are the ones used bv lPun et alJ l|2002l) . 
which are the C NO enriched abundanc e s for H e, C, N and 
O, taken from "Lundavist & Fransson (1996^), and LMC 
abundances from iRussell fc Dopita (^1992. ') for the other 



elemen ts, with exception to Si which is from I Welt v et all 

Because the focus in this paper is on the coronal lines 
we do not attempt to calculate the structure of the pho- 
toionization heated gas behind the cooling, shocked gas. 
This region dominates the flux for medium and low ioniza- 
tion lines, like [O III], [Ne III] and lower ionization stages. 
The structure of this region, as well as the pre-ionized 
emission lines, will be discussed in a separate paper. Here 
we only note that the low and intermediate ionization 
lines, like [O III], [Ne V], [Ar V], and [Fe VII], all have 
narrow components at zero velocity. The higher ioniza- 
tion lines, [Fe X], [Fe XI], and [Fe XIV], however, lack 
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Fig. 7. Evolution of fluxes of the high ionization hues diative shocks as function of the shock velocity . The fluxes 



compared to that of the X-ray flux from .Park ct al. t2005) . 



are given relative to the energy flux, Fshocit = PKhock/^j 
through the shock. 



such a component. This indicates that pre- ionization of 
the unshocked gas does not reach more than the former 
stages . Models similar to those in lLundavist fc FranssoiJ 
l)l99(i() show that gas ionized by the initial UV/soft X- 
ray flash at shock break out with densities ^ 10^ cm~'^ 
is at these epochs mainly found in Fe 1-11, with ^ 10"'^ 
in Fe III. Only in gas with densities ^ 10'^ cm^"^ can ions 
like Fe VII be found. 

The collision strengths of the impor tant coronal lines 
are mainly from t he Chianti data base ijPere et al.lll997( 
iLandi et alJl2006^ . but in some cases supplemented with 
more recent calculations. The critical densities for coUi- 
sional de-excitation of the [Fe X-XIV] lines are in the range 
10^ — 10^ cm""^, and coUisional destruction is therefore not 
important for these. Consequently, the relative line ratios 
are not sensitive to the exact value of the density. Here we 
take a value of ngpot = 10'' cm~^ for the pre-shock density. 

In Fig. |S1 we show the strengths of the most impor- 
tant optical, high ionization lines as a function of the 
shock velocity relative to the kinetic flux through the 
shock, pV^/2. From this we see that at shock velocities 
^ 400 km s^' the line ratios vary strongly with velocity, 
while above this velocity they are, with the exceptions of 
[S XII] A7611.0 and [Ar XIV] A4412.2, nearly constant. 
This is explained by the fact that above ~ 400 km s^' the 
most abundant ions in the immediate post-shock gas are 
ions more highly ionized than Ne V and Fe XIV. For a ra- 
diative shock with velocity higher than ~ 400 km s^^ the 
fraction of the cooling region above '--^ 2 x 10® K, which 
corresponds to ionization stages higher than these ions, 
is largely irrelevant to the flux of the coronal lines, and 
the line ratios will therefore be insensitive to the shock 
velocity. 

In Fig. El we show the ionization and temperature 
structure of iron of a radiative shock with velocity 
350 km s^^. Immediately behind the shock in a region 
too small to be well resolved on this scale, ~ 5 x 10^'^ cm, 



the ions adjust from their pre-ionization values to their 
equilibrium abundances. Behind this region Fe IX- XIV all 
have high abundances in most of the shocked gas. As the 
gas cools the temperature and ionization slowly decreases. 
When the temperature has fallen below ~ 10® K cooling 
increases catastrophically and drops to (1—3) x 10'' K. This 
region is too thin to be resolved in the plot, but is very 
important for the observed emission lines (see below). The 
different lines are therefore direct probes of the tempera- 
ture interval where they have their maximum abundance. 




Distance from shock (10"^ cm) 

Fig. 9. Ionization structure for iron and temperature be- 
hind a radiative shock with velocity 350 km s~' and pre- 
shock density lO'' cni^^^. 

To more clearly illustrate the sensitivity of the dif- 
ferent lines to the temperature we show in Fig. IIUI the 
contribution to the luminosity per logarithmic tempera- 
ture interval of each line as function of temperature. This 
clearly shows how the different lines trace different tem- 
perature ranges behind the shock, with the higher ion- 
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ization stages naturally more sensitive to the higher tem- 
peratures. Consequently, the [Fe XI-XIV] lines arise at 
(1 - 2) X 10^ K, while the [Ne V], [Ar V] and [Fe VII] 
lines arise at 3 x 10^ K in a geometrically narrow, but 
important, region, where the temperature falls from 10^ K 
to 10^ K. Because the whole post-shock region is in near 
pressure equilibrium, the density is ne(T'e) ~ 4nspotTs/re- 
The [Ne V] and [Fe VII] lines e.g., arise at Te - 3 x 10^ 
K, corresponding to ^ 24nspot ~ 2.4 x 10^ cm~^, for a 
pre-shock density of rigpot — 10^ cm^'^. The critical den- 
sities for these lines are 10^ cm~'^, and coUisional de- 
excitation is therefore clearly not important for these. 

For completeness Fig. 1101 also shows the [N II], [O III] 
and [Ne III] lines. In the coUisionally ionized part of the 
shock they each trace different temperature intervals from 
~ 10'' K to 4 X 10^ K. They all have, however, more im- 
portant contributions from the photoionized, cooling gas, 
not included in these models. 

To compare the models with our observations we relate 
all hne fluxes to that of [Fe XIV] A5302.9 at the different 
epochs. These ratios are given for the different epochs in 
Table In the same Table we also give the shock velocity 
corresponding to the different line ratios, taken from Fig. 
ISI As can be seen, the range of shock velocities is for all 
line ratios surprisingly small, 310 — 390 km s~^. The corre- 
sponding temperature behind the shock is (1.6 — 2.5) x 10^ 
K. This comparison assumes that the emitting area in each 
line is the same. If there are shocks with a range of veloc- 
ities this assumption may be questionable. Because of the 
rapid variation of the line ratios with velocity, this is, how- 
ever, not likely to affect the derived velocities by a large 
factor. 

There is some tendency for lines from ions with lower 
ionization potential, e.g., Ne V, Ar V and Fe VII, to give 
a somewhat lower velocity than the lines from high ioniza- 
tion ions like Fe X and Fe XI. The fact that we do not see 
[Ar XIV] A4412.2, at a flux level of ^ 15% of the [Fe XIV] 
line is consistent with a shock velocity ^ 500 km s~^. 
The [Ca XV] A5694.4 line gives a similar upper limit of 
- 500 km s-i (Fig.lHl). 

Finally, we note that although only barely covered by 
these observations, one also expects the [Fe XIII] A3388.5 
and AA10746.8 - 10797.9 lines to be detectable in the UV 
and near-IR, respectively. In particular, the IR lines are 
expected to have a strength comparable to that of [Fe X] 
A6374.5 (see Fig. EJ. The A3388.5 line is within our ob- 
served range, but the spectrum is here noisy due to the 
decreasing sensitivity of the instrument and the influence 
of the atmosphere. We estimate a maximum strength com- 
parable to that of [Fe XIV] A5302.9. Because the expected 
flux is ~ 10% of this line, this limit is not very useful. The 
IR lines are outside the range. 

5. Discussion 

Although we get good agreement for a rather narrow 
range of shock velocities, it is clear that the dynamics 
is more complex. This can be seen in the hydrodynamic 



[FeXIV] 5302 

n=10*cm'^ 




Temperature (K) 

Fig. 10. Contribution to the luminosity per logarithmic 
temperature interval to different lines behind a shock with 
velocity 350 km s~^. 

simulations of the collision performed by several groups 
iILuo et al.ll994HMasai fc Nomotoll994l:lBorkowski et al. I 
Il997l) . In particular. iBorkowski et al. I I^1997^ have mod- 
eled the structure and X-ray emission from the impact 
with 2-D hydro simulations. The interplay of the shock 
propagating into the dense ring gas and reflected shocks 
from the ejecta-CSM gives a time dependence, as well 
as spatial dependence, which is complex and only can be 
described by detailed multi-dimensional hydrodynamical 
simulations. Although the Borkowski et. al. simulations 
used a simplified model for the ring structure, they give a 
good idea of the processes and complications involved in 
modelling the emission. 

A further complication is that the observed spectra 
are a convolution of several hot spots seen with different 
projected velocities. In addition, the number of spots is 
increasing. In these respects our observations have similar 
limitations as the Chandra X-rays observations, although 
our spectral resolution is much better. Spatially resolved 
spectroscopy would be of obvious importance. 

In Sect. 3 we found that the [Fe XIV] A5302.9 line 
extends to ~ 350 km s~^, while the lower ionization lines 
only extend to ^ 250 km s~^. This may indicate that the 
low ionization lines have a dominant contribution from 
tangential shocks with lower velocity, while the high ion- 
ization lines also come from shocks perpendicular to the 
bla st wave, havin g a high velocity. 

iMichael et al.l 112002') find that the relative fluxes of the 
lines observed with the grating spectrometer on Chandra 
in October 1999 are best fitted with a shock temperature 
of kTe - 2.9 ±0.4 keV. The lower dispersion CCD spectra 
obtained in December 2000 indicated a marginally lower 
temperature, kT^ ~ 2.7 ± 0.2 keV. In addition, they flnd 
that they get a slightly better fit to the spectrum with 
an additional softer component with fcTe ~ 0.4 ± 0.1 keV. 
The flux contribution from this was ~ 7%. It is worth 
noting that even the high temperature component has an 
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Table 3. Reddening corrected fluxes relative to [Fe XIV] A5303 and corresponding shock velocities ( km s^^). 



Line Epoch 2 Epoch 3 Epoch 4 Epoch 5 



[Fe VII] A6087.0 


_ 


0.19 


0.10 


0.094 






316 


330 


332 


[Fe X] A6374.5 


0.65 


0.68 


0.54 


0.57 




369 


366 


391 


381 


[Fe XI] A7892.0 


1.05 


0.63 


0.68 


0.53 




326 


340 


338 


348 


[Ne V] A3425.9 




1.74 


1.00 


0.83 






310 


318 


321 


[Ar V] A7005.7 




0.076 


0.025 


0.021 






313 


336 


339 



electron temperature that is considerably lower than the 
shock temperature, corresponding to equipartition, ~ 17 

keV. This is in common with several other SNRs. 

Due to the increase in the X-ray flux IZhekov et al.l 

l)200 5. 2006) were able to repeat this study with consider- 
ably higher S/N. Zhekov et. al. find that they get a good 
fit to the X-ray line ratios, as well as the widths, for a 
model with a range of temperature, 0.15 — 4 keV, corre- 
sponding to Vs = 340 — 1 700 km s^^. The lower part of 
this range is similar to that observed in the optical lines, 
an d it is very like l y tha t they arise from the same region. 
As IZhekov et al.l |2005(l remark, the upper range of ve- 
locities is, however, likely to come from adiabatic shocks, 
probably reflected from the impact on the protrusions. If 
they are adiabatic and with velocity ^ 400 km s^^, corre- 
sponding to a temperature of ^ 3 x 10^ K, they will not 
produce optical emission. 

We propose that the regions that produce the soft X- 
ray flux also account for most of the emission in the [Fe X] , 
[Fe XI] and [Fc XIV] lines. This is supported by the sim- 
ilar evolution of the flux between these two components 
(Fig. [71), and, as mentioned above, the comparable shock 
velocities required to produce them, ^ 350 km s~^ (see 
also Table 3). However, in order to explain the slower evo- 
lution of the intermediate ionization lines such as [Ne V] , 
[Ar V] and [Fe VII], we do require that some fraction of 
the [Fe X-XIV] fluxes arise from adiabatic shocks which 
would not be contributing to the intermediate and low ion- 
ization lines. The sensitivity of the cooling time scales to 
the shock velocities and the density of the hot spots (Eq. 
121) makes this a plausible scenario. The somewhat narrower 
line widths that the intermediate ionization lines exhibit 
also support this argumentation. 

6. Conclusions and summary 

The most important result of this paper is the discovery 
of a number of high ionization lines coming from gas of 



temperatures up to ~ 2 x 10^ K. The rapidly increasing 
flux of these is well correlated to the flux in the soft X- 
rays, and offer a complimentary view of the interaction of 
the ejecta and the ring material. The large range in ion- 
ization stage with similar line profiles shows that most of 
the emission is coming from radiative shocks with veloc- 
ity 310 — 390 km s~^, although a fraction of the coronal 
line emission may originate in adiabatic shocks. The shock 
velocity we find from our spectral modeling is consistent 
with the width of the lines. We hope in the future to do a 
more detailed modeling of the spectra, including both the 
lower ionization lines from the photoionized, cooling gas, 
as well as the X-ray lines. Also a more detailed modeling 
of the hydrodynamics o f the explosion along the lines of 
iBorkowski et al. 1 ljl997l) with these extra constraints from 
the line strengths, as well as the line profiles would be 
highly interesting. Continued monitoring of the ring colli- 
sion is of obvious high importance. The increasing flux, as 
well as new instruments using adaptive optics, will enable 
us to study this unique collision in even more detail. 
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